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ABSTRACT 

The  probability  distribution  of  the  minimum  of  the  dis- 
tances from  a  randomly  occurring  trouble  spot  to  n  carriers 
patrolling  along  a  line  of  length  L  is  analyzed.   Two  approxi- 
mations, a  Poisson  process  and  a  fixed  lattice  of  equally 
spaced  points,  which  bound  the  analytic  model  are  developed 
and  their  usefulness  and  limitations  are  discussed.   It  is 
hypothesized  that  a  two  dimensional  Poisson  field  and  a  two 
dimensional  lattice  of  fixed  points  will  form  upper  and  lower 
bounds  for  the  more  algebraically  tedious   two-dimensional 
area  patrol  model.   The  hypothetical  bounding  distributions 
are  developed  for  n  units  patrolling  an  area.   Finally  a 
one  dimensional  radius  of  influence  model  is  developed  which 
quantifies  the  contribution  that  the  effective  operational 
radius  of  a  carrier  aiirwing  makes  to  the  initial  minimum 
distance  analytic  model. 
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I.   INTRODUCTION 

The  Commander  in  Chief  Pacific  has  the  basic  responsi- 
bility of  exercising  control  of  the  military  forces  of  the 
United  States  from  the  west  coast  of  North  America  to  the 
east  coasts  of  Asia  and  Africa.   Since  most  of  this  vast 
area  is  either  ocean  or  accessible  by  sea  lines  of  communi- 
cation, the  U.S.  Navy  has  a  critical  role  in  the  maintenance 
of  a  U.S.  military  presence  in  the  Pacific  Command  area. 

The  problem  addressed  in  this  thesis  was  posed  by  Vice 
Admiral  E.P.  Aurand,  USN  to  Karl  Eulenstein^  Senior  Analyst  in 
Section  J  77  at  CINCPAC.   As  COMTKIRDFLT,  Admiral  Aurand 
wanted  to  know  how  many  aircraft  carriers  would  be  required 
to  patrol  the  western  littoral  of  the  Pacific  Ocean  Basen 
as  well  as  the  littorals  of  the  Indian  Ocean  Basen  in  order 
that  at  least  one  patrolling  carrier  would  be  within  a 
command  specified  distance  of  a  randomly  occurring  trouble 
spot  with  at  least  a  command  specified  probability.   More 
foirmally,  what  is  the  probability  distribution  of  the  minimum 
of  the  distances  between  n  carriers  and  a  single  random 
trouble  spot  somewhere  along  the  patrolled  coast 

This  thesis  addresses  the  analytical  expression  of. 


P^[D  =  min(D^,D2,.  .  .,D^)  £d]  (1) 
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where: 

n   is  the  number  of  carriers  on  station 

D, ,D^  ...  D   are  the  distances  from  each  carrier  to 
12       n 

a  random  trouble  spot 

d   is  a  command  specified  distance  measured  from, 
the  trouble  spot. 

The  model  allows  the  determination  of  the  minimum  n  required 
in  expression  (1)  such  that  the  probability  is  at  least  as 
large  as  p.   Other  related  inferences  and  applications  are 
discussed. 

In  the  following  a  one-dimensional  model  and  two  approxi- 
mations that  give  useful  bounds  to  it  are  developed.   Then 
the  analogous  two-dimensional  bounding  models  are  considered. 
Finally  the  concept  of  a  "sphere  of  influence"  is  incorporated 
into  the  one-dimensional  model. 
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II.   INITIAL  ANALYTIC  SOLUTION 

A.  ASSUMPTIONS 

To  increase  the  tractability  of  the  analysis  several 
simplifying  assumptions  were  made: 

(1)  The  complex  curvilinear  coastlines  v/hich  mark  the 
western  boundry  of  CINCPAC'S  area  of  responsibility  are 
approximated  by  a  straightline  of  length  L. 

(2)  Each  of  the  n  carriers  were  assigned  individual 
patrol  line  segments  of  length  —  miles. 

(3)  A  carrier's  position  is  assumed  to  be  uniformly 
distributed  on  its  assigned  patrol  line  segment. 

(4)  A  trouble  spot  will  appear  at  random  along  the 
coast  line. 

(5)  An  individual  carrier  is  considered  as  a  point. 
(The  effects  the  carrier  airwing ' s  radius  of  influence  are 
analysed  in  Chapter  V.) 

(6)  Each  carrier's  location  within  its  patrol  segment 
is  independent  of  the  locations  of  the  other  n-1  carriers. 
The  location  of  a  trouble  spot  T  is  independent  of  the  loca- 
tions of  the  carriers.  See  Figure  1. 

B .  DEVELOPMENT 

The  Mean  Value  Theorem  for  Integrals  is  utilized  in  the 
derivation  of  the  analytical  results  detailed  below.   This 
theorem  shows  one  may  approximate  a  probability  density 
function  at  a  point  h  as 
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Schematic  and  Symbology  for  the  Analytic  Model 


D.   

w —  Da  -> 


Xa 


<-^3  -H 


D^ 


^.-1 


t — i: 


X3  Xti-I  X-n 

-^M 1  '       I   I  .    If 


X. ,    i  =  1,2,  ...,  n    the  location  of  carrier  i  within 

line  segment  i 

D.,  i  =  1,2,  ...,  n    the  distance  to  carrier  i  measured 

from  a  randomly  occurring  trouble 
spot  t 

L  length  of  line  to  be  patrolled 

n  the  number  of  carriers  assigned  to 

patrol  a  line  of  length  L 

t  randomly  occurring  trouble  spot 


Notationally  this  paper  will  follow  the  custom  of  using 
uppercase  letters  for  random  variables  and  lowercase  letters 
for  parameters  and  outcomes  with  the  exception  of  the 
parameter  L  which  will  denote  the  length  of  the  line  to 
be  patrolled. 


Figure  1 
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f(h)Ah   =   lim   P [X e  (h,h+Ah) ] 
Ah-^0  + 


Examination  of  the  geometry  of  the  problem  suggests 
that  the  development  of  the  cumulative  distribution  function, 
F_(d)    =   P  [D  =  min  (D,  ,  .  .  . /D  )  <_d]   be  segmented  as  follows 

(1)  Left  end  segment  such  that  T  e  (0'  ^  ^ 

(2)  Middle  segments  such  that  T  z  {- — / — )  for  some 

i  =  2 , 3/  . - . ,  n-1 

(3)  Right  end  segment  such  that  T  e   (  ^^"''•^^,L) 


Xaj^JJ'^^-"-!  >•  ,  _1 


O  -^  (r<-l)L  i_ 

n  -n 


^sy-         -   —  -^^,„ ^^ 


Figure  2.   Line  segmentation 

The  general  analytic  result  was  developed  by  conditioning 
on  the  segment  containing  the  target.   As  a  result  of  the 
assumption  that  the  trouble  spot  will  appear  at  random  along 
the  line  (0,L),  the  target  is  located  in  segment  i  with 


The  cumulative  distribution  function  which  will  be 
developed  is  actually  a  function  of  the  minimum  of  the 
distances,   d;  the  length  of  the  patrolled  line,  L;  and  the 
number  of  patrolling  units,  n; Fq (d; L,n) .   For  simplicity 
this  terminology  will  be  abbreviated  as  F_(d). 
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probability  1/n;  that  is. 


P[T£  [^^"j;^^.^]]   =  i-;   for   i  =  1,2,  ...,  n. 

(2) 
Given  that  T  e  [- =-^ — /~7-]  /  we  further  condition  on  T  being 

■Li       Li 

located  in  the  first  or  second  half  of  this  line  segment. 
Let  D..^  =  |T-X._^|,  D.  =  |T-X.|,  D.^^  =  I T  -  X^^.^  1  ,  and 
D  =  [min (D. _,  ,D .  D .  , ) ]  ;  where  D  is  the  random  variable 
whose  cdf  is  being  developed.   The  conditional  pdf  of  D 
given  T  is  unconditional  over  T  (being  located  in  a  specific 
half  of  a  line  segment)  and  integrated  over  the  range  (0,d) 
to  determine  the  conditioned  cumulative  distribution  function 
(cdf)  of  D  =  [min  (D .  _-,  ,D  ■  ,D  .  -  )  ]  given  that  T  is  an  element 
of  a  specific  line  segment  i. 
1.   Left-hand  Segment 

The  assumptions  made  in  developing  this  portion  of 
the  model  are : 

(1)  T  ~ Uniform  (0,L),  so  T  s  (0,-)  with  probability 

n 

1/n 

(2)  X^  -  Uniform  (0/^  ) 

(3)  X^  ~  Uniform  (  ^  ,  i^  ) 

^  n   n 

(4)  X-,,  X^  /  T  are  independent  random  variables, 
a.   Conditional  Probability  Density  Function 

Let  S(i)  be  the  event  that  the  random  variable 

T  occurs  in  the  i    segment  {- — , — )  of  the  line;  then 

n    n 

P[S(i) ]  =  1/n. 
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Consider  the  case  t  £  ^^ 'Tn^  '    ^^ 


D   =   D^   =   |T-Xj_|, 


and  the  maximuin  value  of  D  is  ^  ~  ^*  ^^-^  0  <_  d  <_  t. 


fj^j^^g(^j  (d|t)Ad   =   P[D  £  (d,d+Ad)] 


=   P[X,  £  (t-d,t-d-Ad) or  X^  £  (t+d, t+d+Ad) ] 


=   2Adn/L  . 


For  t  <_  d  <_  —  -  t. 


P[D£  (d,d+Ad)]   =   P  [X^  £  (t+d,  t+d+Ad)  ]   =   Adn/L. 


Consider  next  the  case  t  £  (t^— ,— )  ,  so  D  =  min(D,  ,D^)  and 

2n  n  '  1 '  2 

D  £  t.   For  0  <_  d  <_ t. 


P[D  £  (d,d+Ad)]   =   P[X^  £  (t-d-Ad,t-d) 


or   X-^  £  (t+d,  t+d+Ad)  ]   =   2Adn/L 


16 


For  ^  -t  1  d  £  t. 


P[D  £  (d,d+Ad)]   =   P[{X^  £  (t-d-Ad,t-d)  and  Y.^    >   t+d+Ad} 


or  {X^    <   t-d-Ad  and  X^    e    (t+d, t+d+Ad) } 


(Adn/L)  .  (^-  t  -  d)  .^+  (  (t-d)  .^)  .  (^: 


Ad  (?•)  ^  [—  -  t  -  d  -  Ad  +  (t-d-Ad)  ] 
J-i    n 


2Adn/L  .  [1  -  alil^l  . 


Thus  it  is  seen  that: 


2n/L  ;     0  1  ^  1  t 

fDlT^^^l^)   =  for  t  £  (0,^) 

n/L   ;     t  <_  d  £  (|)  -  t 


^  ;  0  <  d  <  C^)  -  t 

jj  n 


f  (l-f);   (^)-tidit 


for   t  £  (oTTf-) 
2n  n 


(3) 
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b.   Unconditioning  Over  te  ^'^'Tn^    ^^^ 

t  £  (^/-)  Given  S(l) 
2n  n 

The  conditional  probability  density  function  is 
graphed  below,  where  values  of  f_^|n,(d|t)  are  shown  for  various 
regions  of  the  d-t  plane. 


Figure  3.   Unconditional  probability  density  function 
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We  now  find  the  conditional  distribution  of  D, 
given  S  (.1)  ,  by  taking  the  expectation  of  f  t-,  i  m  gn)  ^^1^^ 
with  respect  to  the  conditional  distribution  of  T  given 
S(l).   Thus  unconditioning  over  T  (given  S(l))  has  the 
mathematical  form: 

L/n 

^D|S(1)^-^^   =    /    ^DlT,S(l)(^lt)  fT|S(l)(t^  ^t  • 

0 

For  0  <  d  <  L/2n, 


i-d 


n   n  ^,  ,   ^      2n   n 

I 

0 


fDis,i,<^>    =     /    i-r^t  W        X-  E^^ 


L 

n 
.       r         2n   -,   nds    n  •,, 

n 


=   (f)2d  +  2(f)2(^-2d)  +2(^)^l-^)(^-^+d) 
Li  l    n         L       L   n   n 


=   (r)^[—  -  d  -  2nd^/L]  . 
ij    n 


For  ^  <  d  <  -  , 
2n  —   —  n 


L   ,  L 

d  — 

n  n 


fD|S(l)(^)   =    /      L'L^^^   /    — (^-— ^'l^^ 

0  d 


=   (g)^  (E-^)^2(g)2(l-i^)(^-d) 


=   (r)  ^  (—  -  5d  +  2d^n/L)  . 

Ij     n 

19 


Therefore  in  the  lefthand  segment. 


fD|S(l)(^)   = 


(r)^[— -d-  2nd^/L]  ;    for   0  <  d  <  ^ 
ii    n  —   —  Zn 


(^)^[^-5d  +  2d^n/L];   for   ^  <  d  <  ^ 
ii    n  2n  —   —  n 


(4) 


c.   The  Conditional  Cumulative  Distribution  Function 
The  cumulative  distribution  function  (cdf)  of 
D  conditioned  by  S(l)  in  accordance  with  expression  (2)  is 
calculated  by  integrating  the  pdf  over  the  range  (0,d)  of 
D: 

d 


For  0  <  d  <  7^  , 
—   —  zn 


F^|3(,)(d)   =    /  (^)2t|i-h-2nhVL]  dh 

0 


2       3 

,n>  2  r2Ld    d     2nd  ,      , 


D  T^2n^      24 
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For   ^  <    d    <   ^  , 
2n   —        —  n 


d 


L/2n 


,n.  2,3Ld      5d^  .   2nd^       (Ili\  (J:L\  ^  fl\  i  ^  \^ 
^L^     ^    n     "     2  3L     "  ^  n  ^  4n^        ^2^  ^2n^ 

^3L^  4n^  24 


(IL)  2  rlM  _  1^  +  2nd"^  _  23L^      19_, 

^L^     ^    n  2  3L        o.    2  ^  24^     * 

24n 


For  the  left-hand  segment. 


2      3 
-nx2r2Ld   d    2nd  ,  j=    a  ^  j  ^  L 

L     n    2    3L  —   —  2n 

2      3      2 
„       ,..  19  ,  ,n.2r3Ld  -  5d  ,  2nd    23L  -,   ^    L  ^  j  ^  L 

^D|S(1)^^)   =     2A^^L^     ^-R 2-^-31- -:;7-T^  '  f°^  2E-^  ^E 

'  z4n 

1  ;  for  d  >  -  . 


n 


(5) 


2 .   Middle  Line  Segments 

We  next  develop  the  cumulative  distribution  functions 
for  the  minimum  of  the  distances  between  a  random  trouble  spot 
and  n-2  carriers  assigned  to  patrol  stations  i  =  2,3,  ...,  n-1; 
stations  which  have  at  least  one  other  patrol  unit  on  either 
side  of  their  patrol  segment. 
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The  middle  line  segment  conditional  cdf  is  identical 
for  each  of  the  "middle"  patrol  stations  and  is  distinguished 
from  that  associated  with  either  of  the  end  patrol  stations 
as  follows: 

(1)   end  segments 


D  =  min(D, /D„)   or   D  =  min(D   , ,D  ) 
i   2  n— i   n 


(2)   middle  line  segments 


D  =  min (D. _, ,D . ,D.  , )   for   i  =  2,3,  ...,  n-1  . 

L  2L 
A  typical  middle  line  segment  is  segment  2,  (  —  , — )  . 

The  assumptions  for  developing  this  portion  of  the  model 

are: 

(1)  T  -  Uniform  (0,L),  so  T  £(-,—)  with  probability 

n   n 

1 
n' 

(2)  X,  -  Uniform  (0,-)  , 

1  n 

(3)  X^  -  Uniform  (^  ,%  , 

z  n   n 

(4)  X^  -  Uniform  (%/^) 

J  n   n 

(5)  X,  ,X„,X-.,  and  T  are  independent  random  variables 
a.   Conditional  Probability  Density  Function 


^    'J  X.   /L  3L, 

Consider  a  given  t  z  ^n'7~'  '    ^° 


°i  "   |t-x^|  , 


D2   =   |t-X2|, 


D    =   min{D,,D2}  / 
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2L  L 

and   the  maximiim  value   of   D   is 1.     For    0    <    d   <    t  -    —  / 

n  —        —  n 


^DlT,S(2)  ^^1^^^^      ~      P[D  £   (d,d+Ad)] 


=      P[X2    £     (t-d,t-d-Ad)    or   X^    z    (t+d, t+d+Ad) ]       =      2Adn/L, 


For    t   -   -   <    d    <    ^  -    t    , 
n  —        —     n 


P [D  £  (d,d+Ad) J       =      P[{X^e   (t-d,t-d-Ad)    and   X^>  (t+d+Ad) } 


or  {X2  £  (t+d, t+d+Ad)  and  X^  <  (t-d-Ad)}] 


dn    ^2L  _  (t+d+Ad)]  •  ?  +  ^  [t-d-Ad]-? 


L      n  ''L     L'       'L 


2Adn^  .L 


—  [^  -   d   -   Ad] 


^    [1  -  (d+Ad)^] 


3L  2L 
Consider  next  a  given  te  [t-i — )/  so  D^  =  !  t  -  X^  I  ,  D-,  =  1 1  -  X-,  L 

^  2nn  2  '  2'3         '  3' 

D   =   min(D2/D^)    and   the   maximuin  value   of    D    is    t  -— • 
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2L 
For    0    <    d    <    —  -    t    , 
—        —     n 


P[D    £(d,d+Ad)]       =      PiX^e  (t+d,t+d+Ad)    or    X^  e   (t-d, t-d-Ad) ] 


2Adn 
L     • 


For   ^  -t    <    d    <    t  -^, 

n  —       —  n' 


P[D£   (d,d+Ad)]       =      PiiX^e  (t-d, t-d-Ad)    and   X^   >  (t+d    Ad)}    or 


{X^  £  (t+d,t+d+Ad)  and   X^  <  (t-d  -Ad)  }] 


Adn   r3L   /a.^J_L^  j\  iH  ^  Adn  , ,  ,  .,  Lxn 

— p—  •  [-— -  (t+d+Ad)h=-+-^— (t-d-Ad--):^ 

L     n  L    L         n  L 


2Adn  rL    J    .J, 
j-[-   -  d  -  Ad] 


2  Adn 


[1  -  (d+Ad)n/L] . 


Thus  it  is  seen  that: 


^d|t,S(2)  ^^^^^ 


2n 
L  ' 


2n  r,         dn. 


2n 
L 


2n  r,    nd-, 
L  ^-^  "  L  ^  ' 


0  <  d  1  t-  - 


^   L   ^    2L   , 
t--<d<  — --t 

n  —  —  n 


0<d£  ^-t 


2i-  t<d<t-^ 

n     —       n 


^^  ^n  '2^) 


,  ,  .3L  2L. 
2n  n 


(6) 
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L  3L 
b.   Unconditioning  Over  t  e  {—f-^r-)    And 

^T  2T  ^  ^^ 

t  £  (|^'^)  Given  S(2) 
zn  .   n 

The  conditional  probability  density  function  is 

graphed  below,  where  values  of  f_^|rn(dlt)  are  given  for 

various  regions  of  the  d-t  plane. 


t  -^^ 


Figure  4.   Unconditional  probability  density  function 


The  conditional  distribution  of  D,  given  S(2), 
is  now  found  by  taking  the  expectation  of  f t-,  i  m  c/?)  ^^^^^  with 
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respect  to  the  conditional  distribution  of  T  given,  S(2), 
Thus,  unconditioning  over  T  (given  S(2))  has  the  following 
mathematical  form: 


2L/n 

^D|S(2)^^^   =     /     ^D|T,S(2) (^1^^  ^T|S(2)^^^  ^^ * 

L 
n 

For  0  <  d  <  T—  , 
—   —  2n 

d+(L/n)  ^d 

L/n  d+- 

n 


2L/n 


f  2n  ,  -   nd^  /n,  , , 

/    -r^i -T-)  ^r^^t 


(2L/n)-d 


n«2,,   nd,  ,,  ,  L   L>  .  ^,n^2,2L         L 
^)  (i._)  (d  +  ---)  +  2(j-)  (_-d-d-- 


+  2(^)^1-^)  (^  -  ^  +  d) 
L       L    n    n 


2  (2.)  2  (L  _  2nd^ 


For  ^  <  d  <  -, 
2n  —   —  n' 


(2L.3,  2L 

"     2n,,   nd,  ,n,  ,,  ,       r  "  2L 


^DlS(2)'^'   =    /        X<l-f"E'^t.       /    ^(l-nd)(^)at 


n 


d.(t, 
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f  l^,,^(d)   =   2(^)2(1. !ldj(^.d-^) 
D I  S  (2)  L       L    n      n 

■    +  2(2.)2(i.£^)(2ii.d-^) 
h  L   n      n 


L    n        L 


Therefore  in  any  middle  line  segment, 


fD|S(2)'^'   = 


2<E''(H-^)^       O^d^^ 


2 

^,nx2  2L    .,   2nd  .     L  ^  -,  ^  L 

2(_)  (_  -4d--^);   2E^  ^  lE  •    (7) 


c.   The  Conditional  Cumulative  Distribution  Function 
The  cdf  of  D  conditioned  by  S(2)  is  calculated 
by  integrating  the  pdf  over  the  range  (0,d)  of  D; 

d 


For  0  <  d  <  T?^  , 

—   —  2n 


0 


3 

_   ^  .n>  2  Ld    2nd  . 
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and 


^D  S(2)  ^2n^      6 


For  2^  <  d  <  -, 


F,|3(,^(d)   =   |.    /2(a)2(2^.4h.^)dh 

Xj 

IE 

3 
_   5    ^  .n.  2  r2Ld  ^.2      2nd    .2L«,L.  ,',/L.2 

^3L^  ^2n^  •■ 

3     2 
5  ,  ^,nv2r2Ld   ^.,2  ,  2nd    7L   , 
=   6  -^  2(^)  [-^-2d   +  -^^-— -^]. 

12n 


For  the  middle  line  segment, 


3 

^  .n^  2  ,Ld   2nd  ^  ^    n  ^  j  ^  L 

3     2 

r,       /jx        5  ,  -,n,2,2Ld   -,,2  ,  2nd    7L   >   ^    L  ^  ,  ^L 
^D|S(2)^^)   =     6^2(^)  (-^-2d  4--3^-__^);  for  27^£dl- 
'  l2n 

1;  for  d  >  -  . 


n 


(8) 


3 .   Right-hand  Segment 

Development  of  the  conditional  cumulative  distribution 
function,  of  the  minimum  distance  for  the  right-hand  segment  is 
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identical  in  methodology  and  form  to  the  results  obtained 
in  analyzing  the  left-hand  segment  in  paragraph  II.B.l  above 
We  therefore  omit  the  details, 

C.   UNCONDITIONAL  CDF  FOR  n  SHIPS  ON  A  LINE  OF  LENGTH  L 

The  cumulative  distribution  functions  derived  thus  far 
are  conditioned  on  S(i)/  i  =  1,2,  ...,  n.   Since  T  ~ Uniform 
(0,L)  and  each  patrol  segment  is  identically  —  units  long, 
the  probability  of  T  occurring  in  any  specific  segment  is , 


P[S.]   =   P[TE:(i^,ii±i^)]   =   [(-ii±lL^-i^]/L   =   i 
1  n    n  n     n^       n 


Therefore  the  unconditioned  cdf  is  derived  from  (5)  and  (8), 
as  follows: 


n 


^D^^)   =   E  I    ^D|S(i)(^)- 


i=l 


For  0  <  d  <  ^  , 
—   —  zn 

2      3  3 

T-i  /  j\  /l\r->r/n\2, 2Ld   d    2nd  >,  ,  ,   ^,r/riN2,_.  ,Ld    2nd  ^  ■,  ■, 

Fj^(d)   =   (k)^;2[{^)  (— --_--^^)]  +  (n-2)[(-)  (2)  (—  -^L-)]} 

2      2  3      3  3 
2nd   nd   _^  4n^d   -  4n-^d" 

^^  "  nL^        3nL^ 

2  2      2  3      3  3 

6L  nd  -  3Lnd   +  4n  d   -  4n-^d-^ 

3L^ 
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and 


p  /Ji_^   -   5  _  _1_    ^   lOn  -  1 
^D^2n^      6    12  "        12n   * 


For  ^  <  d  <  ^, 
zn  —   —  n 

^  ,,.      (lQn-1)  .  ,l^  ror /n,2,3Ld   5d^   2nd^   23L^, 

^D^^)   =     12n   ^^n^^^^^L^     ("H 2 SlT  "  TTT^ 

24n 

3     2 

.  -   ^>  ,^>  ,n«  2  r2Ld   ^  ,2   2nd    7L   ,  -, 

+  (n-2)  (2)  (-)  [-- 2d  +—^7 y]  ) 

L    n  3L    , T  z 

12n 

(lOn-1)  ,  4nd-2d  ,  3nd^-4n^d^  ,  4n'^d^-4n^d^  ,   5    7 
—  + _ -). _ -I- _ 1-     — 


12n        L        ^2  ^^3       12n   6 


22  2  22  32  43         333         3 

12n   L   d-6L   nd+9n   d   L-12n   d   L+4n   d   -4n   d    +L   -nL 


3L   n 


Therefore  the  cumulative  distribution  function  of  the 
minimum  of  the  distances  from  T  to  n  patrol  units  is: 


2  2         2    3         3    3 

6L   nd-3Lnd   +4n   d   -4n   d                ^           ^    ^     ,  ^     L 
^ ;         for      0    <    d  <    TT— 

3L^  -  -  2^ 


F^(d)       = 


22  2  22  32  43  333         3 

12n    L   d-6L   nd+9n   d    L-12n   d   L+4n   d    -4n   d-^+L   -nL-^ 


3L^n 


for  ^  <    d    <   -    . 

2n  —        —  n 


(12) 
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III.   APPROXIMATIONS  TO  ANALYTICAL  RESULTS 

A.  RATIONALE  FOR  APPROXIMATIONS 

The  existance  of  a  simple,  readily  applicable  approxima- 
tion, especially  one  with  well  understood  and  widely  known 
characteristics  of  behavior,  might  be  useful  to  the  military 
commander  and  planner.   This  thesis  investigates  two  such 
approximations  and  evaluates  the  range  of  values  for  which 
each  may  be  an  appropriate  approximation  of  (12). 

B.  POISSON  PROCESS 

The  Poisson  process  as  a  mechanism  for  positioning  n 
units  on  a  line  of  length  L  is  appealing  because  of  the  well- 
developed  and  tractable  body  of  theory  regarding  the  exponen- 
tial distribution  of  the  distances  between  units.   The  limi- 
tations of  the  approximation  will  be  discussed  after  the 
model  has  been  completely  developed. 

1.   General  Properties 

For  the  model  let  n  be  the  average  number  of  ships 
located  on  a  line  of  length  L,  so  A  =  ^  is  the  average  den- 
sity of  ships  per  unit  length  of  L.   Let  W.  be  the  distance 
between  unit  i  and  i-f-1,  for  all  i:   i=l,2,  ...,  n-1.   Then 
for  a  general  Poisson  process  one  has: 


P[N  =  n]   = 


, _ , n   -AL 
AL)   e 

n! 
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The  interarrival  "distances"  (W.)  between  adjacent  ships  are 
independent  exponentially  distributed  random  variables  such 
that 

W.   ~   Exp(l/A) . 

In  what  follows  only  those  cases  where  T  occurs  between 
two  units  will  be  considered  (the  interval  (0,X,)  will  be 
discussed  later) .   Let  the  trouble  spot  T  occur  at  random 
on  the  line  of  length  L  as  illustrated  below. 

K W,  ^ 


V' ^ L_ 1 ^ 

K—  M'  — i^ A12 


Let  Ml  represent  the  distance  from  ^.  to  t,  M2  the  distance 

from  t  to  X.  ,,  then  it  is  easily  seen  that  Ml + M2  =  W. 
1+1  ^  1 

the  distance  from  the  i    unit  to  the  i+1    unit.   The 
memoryless  property  of  the  exponential  distribution  implies 
that  both  Ml  and  M2  are  distributed  as  exponential  (1/X) 
random  variables.   Now  with  D  =  minCMl,M2),  we  have, 


Fj^(d)   =   PlMin(Ml,M2)  £d]   =   PCD<_d). 


=   l-P(D>d)   =   1-P  (Ml  >  d,M2  >  d) 


=   1  -  P  (Ml  >  d)PCM2  >  d). 
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and  since  Ml  and  M2  are  identically  distributed, 


Fj^(d)   =   1  -  [1  -  F^(d)]^ 


Thus  one  has. 


F^(d)   =   1  -  [1  -  (1  -  e  ^'^)]^ 


D 


-2\d 
=   1  -  e      ,   A  >  0  .  (13) 


The  model  developed  above  for  approximating  the 
distribution  of  the  minimum  distance  from  a  random  trouble 
spot  to  the  adjacent  patrol  units  implies  the  cdf . 

2nd 
Fj^(d)   =   1-e    ^;   0£d  (14) 

which  is  an  exponential  2A  cdf . 

2 .   Limitations  To  The  Poisson  Approximation 

The  Poisson  approximation  cannot  guarantee  that  n 
units  will  be  positioned  within  the  (0,L)  limit  of  the  line 
to  be  patrolled.   For  example  there  exists  a  positive  proba- 
bility, however  small,  that  no  units  will  occur  within 
(0,L) . 

The  model  as  developed  does  not  consider  the  specific 
distributions  involved  when  T  occurs  to  the  left  of  the 
location  of  the  first  unit.   The  effect  of  this  simplification 
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is  to  increase  the  size  of  d  for  an  arbitrary  probability 
p  which  results  in  a  weakening  of  the  boundary  condition 
developed  in  III.D. 

The  Poisson  approximation  does  not  preclude  three 
or  more  units  being  very  close  to  each  other. 

C.   FIXED  STATION  MODEL 

A  second  approximation,  the  fixed  station  model/  is  nov/ 
considered.   This  model  assumes  each  patrol  unit  maintains 
a  fixed  station  located  at  the  midpoint  of  the  assigned 
patrol  line  segment  X.  =  (2i-l)L/2n. 


X. 


Aa  A3 

J 1 1— . 


f  2  4- 


1.   Model  Development 

Since  the  distance  between  adjacent  units  is  constant 
(— )   and  the  distance  from  any  unit  to  its  patrol  line  segment 
boundaries  is  L/2n,  the  model  developed  for  an  arbitrary 
patrol  segment  will  be  representative  of  all  segments. 

Consider  segment  1.   Given  S(l),  T  is  distributed 

uniformly  on  (0,—)  and  the  model  being  developed  will  specify 

the  conditional  cdf  of  D  =  minims T,T -•=—}.   That  is, 

2n    '    2n 
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F^l_,,,(d)   =   P(D<d|S(l))   =   P(D  <  dJT  <  ^)P(T  <  ^) 
D|S(i)  —  '  —  '   —  2n     —  2n 

+  P(D  <  d|T  >:^)P(T  >^) 
—  '    zn     —  2n 


p(^.T<dlTl^)i+  P(T-^ldiT>^)i 


P(T  >^-d|T  <^)i  +  P(T  <  d  +  ^|T  >^)i 
2n    '    2n  2       —    2n '    2n  2 


^  ^_L L_ 

_U-p(T<^-diT<2^)]  +  2^ l721^^ ^ 

^-  d 
1[1  -  |iL^]  +  l(2nd. 


2'     L/2n  '    2'  L 


1[1  -  (iiz2^)(^)]  +  i^ 


_1  _  L  -  2nd    nd 
2      2L      L 


^  ;    for  0  £  d  1  ^  .  (15) 


Expression  (15)  is  conditioned  by  the  occurrence  of 

S(l).   However  the  conditional  cdf  given  S(i)  for  each  i; 

i  =  1,2,  ...,  n;  of  D  =  min (D . ,D ., , )  is  identical  to  (15) 

1   1+1 

and  the  summation  of  n  conditional  cdf's  multiplied  by  the 
conditioning  probability  1/n  results  in; 


P(D  <  d)   =   ^  ;    for   0  <  d  <  ^  (16 

—         L  —   —  2n 


which  has  a  median  value  of  -r— .  (17) 

4n 
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2.   Limitations  Of  The  Fixed  Station  Model 

The  obvious  limitation  of  the  model  is  the  loss  of 

realism  in  constraining  a  mobile  unit  to  a  small  geographically 

fixed  station.   As  a  consequence  of  this  constraint  the 

maximum  value  of  D  =  min(D.,D.  ,,)  is  77—  for  the  fixed  station 

1   1+1     2n 

model . 

D.   COMPARISON  OF  ANALYTIC  AND  APPROXIMATE  MODELS 

The  easiest  method  of  contrasting  the  two  approximations 
with  the  original  analytic  model  is  through  a  tabulation  of 
their  individual  cumulative  distribution  functions  for 
specific  values  of  n,  L,  and  D.   In  Table  1  d  is  tabulated 
as  a  percentage  of  the  ratio  — .   This  procedure  allows  one 
to  readily  compare  the  various  models  at  differing  L  and  n 
values . 

For  the  L  =  30000,  n  =  6  case  the  distribution  function 
for  the  analytic  model  is  bounded  below  by  that  of  the  Poisson 
process  and  above  by  that  of  the  fixed  station  model. 

Figure  (5)  is  a  graphic  presentation  of  the  "fit"  of 
the  2  bounding  models  when  compared  to  the  analytic  model 
with  parameters  L  =  30000  and  n  =  6.   The  graph  depicts  the 
residuals  obtained  from  differencing  the  analytic  and 
approximation  models  at  various  values  of  d.   Data  points 
are  from  Table  1. 

For  the  case  under  consideration  the  Poisson  model  median 
value  of  D  (0.3465)  does  not  approximate  the  analytic  model 
median  (0.266)  as  closely  as  does  that  of  the  fixed  station 
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model  (0.25).   Except  for  values  in  the  neighborhood  of 
d  =  .5  the  Poisson  process  model  residuals  are  of  greater 
magnitude  than  are  the  fixed  station  residuals. 

The  Poisson  model  gives  reasonable  (if  somewhat  conserva- 
tive) estimates  of  the  analytic  model  results  within  the 
range  0.0  <_  d  <^  .3. 

For  this  case,  the  fixed  station  approximation  model 
not  only  bounds  the  analytic  model,  from  above,  it  also 
appears  to  be  a  reasonable  approximation  of  the  analytic 
model  over  most  of  the  range  of  d£  (0.0,1.0).   The  fixed 
station  median  value  of  0.25  is  very  close  to  the  analytic 
model  median  value  of  0.266.   Using  a  criterion  that  a 
residual  value  greater  than  0 . 1  indicates  that  the  approxima- 
tion model  is  inappropriate  for  that  value  of  d,  it  is  seen 
that  the  fixed  station  model  is  appropriate  for  all  d  except 
0.4  <  d  <  0.6.   This  result,  while  developed  for  the  case 
n  =  6  and  L  =  30,000  will  hold  in  general  except  in  those 
unusual  cases  where  n  >  L. 

The  results  tabulated  in  Table  1  suggest  that  there 
exists  a  natural  hierarchy  or  ordering  of  the  three  models. 
If  one  fixes  the  P(D<_d)  =  .5  and  then  ranks  the  models 
according  to  the  relative  magnitudes  of  d  required  by  each 
model  to  obtain  this  probability  the  result  is  as  follows: 


Model 

Median 

1. 

Poisson 

0.3465 

2. 

Analytic 

0.266 

3. 

Fixed  station 

0.25 
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Arranged  in  this  manner  the  models  reflect  the  following 
"costs" : 

1.  Increasing  level  of  control  by  higher  operational 
authority 

2.  Stochastically  decreasing  distances  to  random 
targets 

3.  Decreasing  randomness  of  own  units  motion 

4.  In  a  hostile  environment  an  increasing  risk  of  an 
adversary  accurately  determining  the  positions  of 
the  patrol  units 

5.  Decreasing  probability  of  adjacent  units  being  in 
close  proximity  (overlapping  of  coverage)  which  is 
paralleled  by  decreasing  problems  with  mutual  inter- 
ference and  contact  identification. 

The  geometry  of  the  residuals  plotted  in  Figure  5  indi- 
cates that  the  average  of  the  two  boundary  cdf ' s  might  form 
a  more  reasonable  estimate  of  the  analytic  model  than  either 
of  the  boundary  models  does  when  considered  individually. 

The  boundary  model  data  in  Table  1  was  averaged,  sub- 
tracted from  the  analytic  (L  =  30000,  n  =  6)  model  and  the 
residuals  plotted  on  Figure  5  as  the  average  model  residuals 
The  resulting  residuals  are  obviously  a  very  reasonable 
conservative  approximation  and  lower  bound  for  the  analytic 
model,  particularly  for  those  ds  (0,.65). 
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IV.   APPROXIMATIONS  FOR  TWO-DIMENSIONAL  PATROL  AREAS 

The  previous  chapters  have  developed  a  one -dimensional 
patrol  model  and  two  approximations  which  form  upper  and 
lower  bounds  for  the  analytic  model.   It  seems  a  reasonable 
supposition  that  when  the  problem  scenario  is  expanded  to 
include  units  patrolling  a  two-dimensional  area  that  a  two- 
dimensional  Poisson  field  would  form  a  lower  bound  of  the 
two-dimensional  analytic  model  and  that  a  two-dimensional 
lattice  of  fixed  stations  could  reasonably  be  postulated  as 
forming  an  upper  bound  on  the  two-dimensional  analytic  model. 

The  algebra  involved  in  developing  the  two-dimensional 
analytic  model  for  n  randomly  patrolling  units  is  quite 
involved  and  time  constraints  precluded  its  development  in 
this  thesis.   However  the  two-dimensional  approximations  are 
tractable  and  are  presented  here. 

A.   TWO-DIMENSIONAL  POISSON  FIELD 

1.   General  Theorv  For  A  Poisson  Field 

Consider  an  array  of  points  distributed  over  a 
plane  P.   For  each  subset  B  of  P  let  N(B)  be  the  number  of 
points  within  B.   The  array  is  said  to  be  distributed  in 
accordance  with  a  Poisson  process  with  density  z,    if  the 
following  assumptions  are  met: 

(a)  the  number  of  points  in  non-overlapping  areas 
are  independent  random  variables 

(b)  for  any  subset  B  of  finite  area,  N(B)  is  Poisson 
distributed  with  mean  L,h ,    where  b  is  the  area 
of  B.  [Ref.  1] 
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The  dimensional  analysis  of  r   is   ^ POJ-J"^ts ^^ 

length  squared 

Thus  if  N(b)  represents  the  number  of  points  in  a  region 
with  area  b: 


/  ,  V  n   ~Cb 
P(N(b)  =  n)  =   ^^'';,  ^ .  (20) 


The  probability  that  there  are  no  points  within  a  radius  r 
of  a  given  location  is  determined  by 


2 

P(R>r)   =   P  (no  points  within  a  region  of  area  iir  ) 


_    -^r'^'^    '  (21) 

=   e 


2.   Poisson  Field  Model 

a.  Assumptions  and  Symbology 

The  following  assumptions  are  made: 

(1)  The  n  units  are  distributed  stochastically 
by  a  Poisson  process  with  rate  z, . 

(2)  The  trouble  spot  T  is  equally  likely  to 
occur  anywhere  within  the  patrolled  area. 

(3)  The  area  to  be  patrolled  is  L  units  in 
length  and  W  units  in  width  (L  >  W)  . 

b.  Development 

In  accordance  with  the  theory  for  a  Poisson 
field,  c  =  y-Tj  and  following  expression  (21),  the  Poisson 
field  model  is. 
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^2 
P(D<d)   =   1-e 


2 
T     -ird  n/LW  -^-^^ 

=   1  -  e  .  (22) 


c.   Limitations  to  the  Poisson  Field  Model 

The  Poisson  field  model  cannot  guarantee  that 
n  units  will  be  positioned  within  the  patrol  area.   Nor  can 
the  model  preclude  five  or  more  units  from  being  assigned 
positions  which  are  in  very  close  proximity  to  each  other. 

B.   FIXED  LATTICE  OF  POINTS 

The  cdf  of  the  model  which  formed  the  upper  bound  of  the 
one-dimensional  analytic  model  was  a  series  of  patrol 
stations  fixed  in  the  middle  of  each  patrol  line  segment. 
A  reasonable  two-dimensional  approximation  would  seem  to  be 
an  array  of  patrol  stations  fixed  in  the  centers  of  their 
respective  patrol  areas.   It  is  hypothesized  that  this  two- 
dimensional  array  is  an  upper  bound  for  the  two-dimensional 
analytic  model. 

1.   Assumptions 

Assume  n  points  are  located  within  an  L  by  W  (L  >  W) 
area  as  shown  in  Figure  6  below. 

Assume  that  the  location  of  T  is  uniformly  distri- 
buted over  the  L  by  W  rectangle.   The  n  ships  are  to  maintain 
station  at  the  center  of  their  respective  L//n  by  W//n 
areas . 
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Figure  6.   Fixed  lattice  of  points 

2 .   Development 

Since  each  of  the  n  ships  is  centered  in  an  iden- 
tically shaped  area  the  model  developed  for  any  one  unit 
will  hold  true  for  all  units.   The  maximum  value  of  D  is 
the  distance  from  the  center  to  the  corner  of  a  unit's 
patrol  area  as  shown  in  Figure  7. 

a.   Conditional  Cumulative  Distribution  Function 
Define  event  S(i)  as  the  event  Te   (region  i) 
with  P[S(i)]  =  1/n,  and 


^D|S(i)  ^"^^   "   P(D  <  d|S(i))  . 


44 


Patrol  Area  X^  of  Fixed  Lattice  Model 
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length  g   =   2(d^-a^)^/2   ^   2  (d^  -  5^)^/2 

4n 

length  h   =   2(d^-c^)^/2   ^   2  (d^  -  ili)  1/2 

4n' 

length  a   =   -^   ,   length  b   =   (^  +  ^)  1/2 
2/K  4n  ^  4n^ 

length  d  =   d,  the  radius  of  the  expanding  circle. 


length  c 


2/n 


3   =   2  cos"^  #, 

d' 


(f)  =   2  cos"-^  I, 
d 


arc  length  s   =   de  ,   arc  length  t   =   d0 

Ag   =   area  of  the  segment  under  s  but  outside  the  rectangle 
"^'  ^  '2     -1,  w  .      2   w2  1/2^  W  , 


^(ds-ga)   =   d^  cos~l(— ^)  -  (d^-  — 

2d/^         4^ 


(■ 


A^   -   area  of  segment  under  t'"but  outside  the  rectangle 
■■  "  2 


=   2^^^~^c)   = 


d^  cos"^(  ^  )  -  -ii-(d2  -Ll)l/2 
2d/n     2/n      "^^ 


Figure  7 
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Refer  to  Figure  (7)  and  consider  the  case  of 
the  inscribed  circle  of  radius  d  such  that,  0  <_  d  <_  W/2/n, 

2 

P(D^d|S(i))   =  M^)i^)      =   ^  . 

n 

Consider  the  case  of  W/2/n  ±  <^  ±  L/2/n  where 
the  expanding  circle  has  intersected  the  longer  sides  of 
the  patrol  rectangle. 


^D|S(i)"^'   =   ^-  '2A3)(5^). 


As  the  circle  radius  continues  to  increase,  consider 
L/2/n   <   d  <_      (^.  |i)l/2  , 


for  which 


n-rrd  ,_,        .     -,^     x  ,      1 


^D|s(i)«i'    =    Tvr  -  '^A^  ^  2A^><o^: 


Thus  the  conditional  cdf  is , 


^2 

for    0    <    d    <   — 


2/n 
,2'        2nA  „  , 

^D|S(i,'^'       =  ^-^  fo.      —    <d<    — 

rnrd^        2nA^        2nA  ^  2  2 

Tw Lir  -  -BT'    f""^  —  1  <i  1  (4H  ■"  IE' 

2v'n 
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b.   Cumulative  Distribution  Function 

To  uncondition  over  the  event  S(i)  one  needs  to 
take  the  following  summation/ 


n 


^D^^^   =  k      I    ^D|S(i)(^) 


i=l 


=   F 


D|SU)  ^"^^  ' 


since  the  n  conditional  cdf's  are  identical.   Therefore  it 
is  seen  that: 


2 

nird  ^    «    -,     W 


LW  ' 


for  0  <  d  < 


2/n 


nTTd^  -  2nA^  ^     W     .     L 


^D^^)   =    LW -'  for  -::-  <  d  < 


2  2  /n        2  /n 

niTd   -  2nA   -  2nA^        ^  ,^2^2,,^ 

s       t^     L^j^.W^L. 1/2 

rn /  fo^  <  d  <  (i —  +  -r— ) 

LW  ^  / —  4n    4n 

2/n 

(23) 


C.   EVALUATION  OF  THE  APPROPRIATENESS  OF  TWO-DIMENSIONAL 
APPROXIMATING  MODELS 

Expressions  (22)  and  (23)  form  hypothetical  bounds  to 
an  undeveloped  area  patrol  model.   If  the  hypothesis  is 
correct  it  would  seem  reasonable  that  the  two-dimensional 
bounding  functions  should  behave  in  a  manner  similar  to  that 
of  the  one-dimensional  bounding  functions.   That  is  one 
would  not  only  expect  the  fixed  array  model  to  yield  sto- 
chastically larger  results  than  the  Poisson  field  model  but 
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one  would  also  expect  to  see  very  close  agreement  between 
the  two  models  for  small  values  of  d.   Table  2  which  tabu- 
lates the  evaluation  of  the  two  models  for  parameters  n  =  25, 
L  =  30000  and  W  =  20000,  indicates  that  for  this  case  the 
two-dimensional  models  are  indeed  well  behaved. 

In  the  one-dimensional  case  it  was  found  that  the  average 
of  the  boundary  models  was  itself  a  very  reasonably  estimator 
of  the  analytic  model.   If  the  hypothesis  regarding  the  two- 
dimensional  boundary  models  is  correct  then  a  similar  averaging 
process  could  be  considered  a  better  approximation  of  the 
analytic  area  patrol  model  than  either  boundary  model  is 
when  considered  separately.   The  average  of  the  two-dimensional 
boundary  models  is  also  tabulated  in  Table  2. 

Model d  values*  (in  units  of  length) 

0    500    1000    1500    2000    2500    3000    3500    3605 

Poisson  field   1    .032   .123    .255    .408    .559    .692    .798    .818 

fixed  point     1    .033   .131    .295    .524    .733    .920    .998    1.0 
array 

average  1         .032        .127         .275  .466  .646  .806  .898  .909 

*    for   parameters    n   =    25,    L   =    30,000,    W  =    20,000. 


Table  2.   Evaluation  of  Two-dimensional 
approximating  models 
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V.   ANALYTICAL  MODEL  WITH  RADIUS  OF  INFLUENCE 

The  initial  analytic  model,  equation  (12),  was  developed 
under  the  assumption  that  an  aircraft  carrier  be  considered 
a  point  on  its  patrol  line  segment.   This  restriction  effec- 
tively eliminated  from  the  previous  analysis  the  known  ability 
of  a  carrier's  airwing  to  exercise  considerable  influence 
over  the  environment  adjacent  to  the  carrier.   We  now  modify 
the  model  to  allow  an  examination  of  the  contribution  to 
PCD<_d)  made  by  an  airwing  controlling  a  radius  r  on  either 
side  of  the  carrier's  position. 

An  additional  benefit  of  relaxing  the  zero  radius  of 
influence  assumption  may  be  a  quantitative  model  of  how  a 
units'  probability  of  being  within  .a  specified  distance  of 
a  random  target  will  vary  as  the  radius  of  influence  is 
changed. 

A.   ASSUMPTIONS 

All  the  assumptions  made  for  the  initial  analytic  model 
remain  in  effect  with  the  exception  of  assumption  II. A. (5) 
which  constrained  a  carrier  to  being  considered  a  point  on  the 
patrol  line  segment.   In  this  modification  each  carrier  will 
be  considered  to  have  a  radius  of  influence,  r,  on  either  side 
of  the  carrier's  position.   If  the  trouble  spot  T  occurs  within 
±  r  of  the  carrier's  position  the  distance  from  the  carrier 
to  the  trouble  spot  is  defined  as  zero.   Distance  to  the 
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trouble  spot,  T,  will  be  measured  to  the  edge  of  a  carrier's 
radius  of  influence  instead  of  to  the  carrier's  actual 
position  as  was  done  in  Chapter  II. 

There  are  no  restrictions  regarding  the  movement  of  a 
carrier  along  its  patrol  line  segment.   This  means  that  the 
second  carrier  would  not  be  restricted  from  periodically 
exerting  influence  over  a  portion  of  the  first  r  miles  of  the 
third  carrier's  line  segment  (see  X^  in  Figure  (8)). 

B.   MODEL  DEVELOPMENT 

1.   Cumulative  Distribution  Function 
In  the  initial  analytical  model, 


D   =   min(D, /D-, . . . ,D  )   where   D £  (0,— ). 


In  the  current  modification 


H   =  min(H-,  ^H^,  .  .  .  ,H  )   where   H  £  (0,  --r] 


and  since  D.  =  H.  +  r  for  i  =  1,2,  ...,  n  then 


H   = 


D-r     for     r  <  D  <  — 

—   —  n 

0       for     0  <  D  <  r. 


Now  for  h  >  0, 
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the  position  of  the  i    carrier 

the  radius  of  influence  on  either  side 
of  a  carrier's  position 

length  of  the  line  to  be  patrolled 

the  number  of  carriers  patrolling 
the  line 

location  of  a  randomly  occurring 
trouble  spot 

the  distance  measured  from  the  random 
trouble  spot  to  the  nearest  extremity 
of  the  i^n  carrier's  radius  of  influe 


uence 


NOTE:    D.  =  H.  +  r 
1     1 


FIGURE  8.   Schematic  and  Symbology  for  the  Radius  of 
Influence  Model 
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Fjj(h)       =      P[H£h]       =      P[H£h|D  ^rJP  (D  ^r] 

+   P  [H  £h|D  <_r]P  (D  £r] 

=     P  [D-r  <_h|D  >_r]P  [D  >_rj    +   P  [0  £  h|  D  <_  r]  P  [D  ^r] 

P  [D-r  £h,D  ^r]P[D  ^r]         P  [0  <_h,D  £  r]  P  [D  <_  r] 
P[D^r]  "^  pTdTTI 

and  since    [0    <   h]    is    certain,    the    above   reduces    to 


P[r£D<_h+r]    +  P[D£r]      =      Fj^(h+r)    -   F^(r)    +   ^^(r) 


=      F^(h+r)  ;       0   £  h   <_  |-  r. 


(24) 


The  radius  of  influence  model  is 


F„(h)   =   F^(h+r)    for    0  <  h  <  ^  -  r       (25) 

rl  U  —    —  n 


which  is  the  initial  analytic  model  [expression  (12) ] 

evaluated  at  D  =  h  +  r. 

2.   Marginal  Rate  of  Return  With  Respect  to  Radius  of 
Influence  Model 

When  the  radius  of  influence  model  is  differentiated 

with  respect  to  r  the  result  is  the  equivalent  to  evaluating 

the  pdf  of  equation  (12)  at  h+r: 


f^(h)   =   fQ(h  +  r)  (26) 


52 


At  r  =  0  this  equation  is  the  true  pdf  associated 
with  expression  (12) ;  at  values  of  r  >  0  this  function  gives 
the  expected  marginal  return  in  terms  of  P(H<_h)  for  a  unit 
change  in  the  radius  of  influence. 

C.   MODEL  EVALUATION 

The  following  plots  were  developed  to  assist  in  visualizing 
two  specific  realizations  of  the  cumulative  distribution 
function  and  marginal  rate  function  for  the  radius  of 
influence  model. 

L n Type  of  plot Figure 

9 
10 
11 
12 

The  parameters  L  =  30,000  miles  and  n  =  6  ships  were 
chosen  as  being  representative  of  the  imaginary  but  roughly 
realistic  situation  described  in  the  Introduction. 

The  plot  in  Figure  11  appears  to  be  quite  linear  for 

all  values  of  r  up  to  and  in  some  cases  beyond  the  r  =  0 . 0 

median  value  of  1330  miles.   It  is  also  of  interest  to  note 

that  for  large  ratios  of  h  ,    that  for  a  given  h  a  constant 

n 

increase  in  r  results  in  a  near  constant  increase  in  P(H  <_h)  . 
This  characteristic  is  not  observed  in  the  L/n  =  1  case. 
By  the  time  the  distance  from  the  target  has  reached 
3500  miles  the  n  =  6,    L  =  30000  cdf  has  exceeded  0.95.   In 
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other  words,  for  all  values  of  r  ^  0  P(H£3500)  >  .95. 
While  this  may  be  a  counter-intuitive  result  it  is  not 
unreasonable  if  one  remembers  that  a  value  of  h  =  3500  miles 
encompasses  a  7000 +  2r  mile  portion  of  the  total  line. 

D.   EMPLOYING  THE   RADIUS  OF  INFLUENCE  MODEL 

We  consider  some  examples  which  illustrate  use  of  the 
model.   The  operational  commander  may  desire  to  know  the 
minimum  number  of  units  necessary  to  patrol  a  30000  nautical 
mile  coast  such  that  the  nearest  unit  has  at  least  a  60% 
chance  of  being  within  a  2-^   day  transit  of  a  random  spot 
along  the  coast.   Further  he  may  desire  to  know  if  there 
is  any  substantial  advantage  to  loading  the  carriers  with 
Type  A  aircraft  (180  nautical  mile  radius)  as  opposed  to 
loading  with  Type  B  (300  nautical  mile  radius).   The  units 
steam  at  25  kts . 

Solve  equation  (25)  for  n  at  h  =  1500,  r  =  180,  L  =  30000 
and  P(H^h)  >_  .  6 .   The  result  is  n  ~  5 .   Equation  (25)  , 
evaluated  at  n  =  6 ,  L  =  30000,  h  =  1500,  r  =  180  results 
in 


F„(h)   =   0.611   >   0.6, 


as  required. 

The  marginal  rate  of  return  function  (26) ,  with  the 
above  conditions  gives 
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f„(h)   =   0.0003, 


which  indicates  one  could  expect  a  0.0003  increase  in 
F„(h)  for  each  unit  increase  in  r.   The  Type  B  aircraft's 
300  nautical  mile  radius  will  increase  F„(h)  by  approximately 

n 

(300  -  180)  (0.0003)  =  0.036  over  the  value  of  0.611.   (The 
actual  increase  by  solving  expression  (25)  for  both  R  =  180 
and  R  =  300  is  0.03959.) 

As  a  result  an  operational  planner  would  inform  the 
commander  that  it  would  require  a  minimum  of  six  carriers 
to  patrol  the  line.   Additionally,  based  on  an  MOE  of 
P(D<_d),  there  is  no  substantial  advantage  to  be  gained  by 
loading  a  particular  type  aircraft  on  the  six  carriers. 
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MflRGINflL    RATES    FOR 
FIXED    VALUES    OF 
RADIUS    OF    INFLUENCE 
L=    30000         N=    6 
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VI.   EXTENSIONS  OF  THE  ANALYTIC  MODELS 

During  the  course  of  the  development  of  this  thesis 
several  potentially  interesting  extensions  were  discussed. 
Some  of  these  ideas  are  presented  below  as  possible  areas 
of  future  investigation, 

A.   MULTI-DIMENSIONAL  PATROL  AREAS 

The  assumption  limiting  a  carrier  to  patrolling  along 
a  line  is  unnecessarily  restrictive.   A  natural  extension 
would  be  to  expand  the  patrol  line  to  a  patrol  area  of 
length  L  and  width  W  such  that  each  carrier  would  be  assigned 
a  L/n  x  W  patrol  area  as  illustrated  below: 


L 


T 
w 


T,    the  trouble  spot,  would  be  restricted  to  either  occurring 
along  the  (_0,L)  line,  or  perhaps  to  occurring  anywhere  within 
the  (L,W)  area. 

An  extension  which  moves  even  farther  from  the  patrol 
line  model,  but  which  naturally  follows  the  previous  exten- 
sion is  to  consider  a  grid  imposed  on  an  open  ocean  area  with 
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one  carrier  assigned  each  rectangular  patrol  area.   T  would 
occur  randomly  anywhere  within  the  patrol  area  boundaries. 
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Figure  13.   n  carriers  randomly  patrolling  a  2  dimensional  area 

The  radius  of  influence  of  each  carrier  could  also  be 
incorporated  into  the  two-dimensional  models . 

B.   GAME  THEORY 

The  stipulation  that  the  location  of  the  trouble  spot 
is  a  random  variable  with  a  particular  probability  distribu- 
tion may  be  realistic  in  most  situations  involving  a  naval 
presence  mission.  Where  an  adversary  relationship  can  be 
postulated  this  stochastic  assumption  as  well  as  the  assump- 
tion of  independence  between  the  location  of  T  and  the 
positions  of  the  carrier  units  (X^'s)  may  be  untenable.   If 
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the  adversary  is  capable  of  some  degree  of  knowledge  of  the 
carrier's  positions,  through  radio  direction  finding, 
shadowing  or  satellite  reconnaissance,  he  is  capable  of 
choosing  T  so  as  to  maximize  his  return  (perhaps  by  maxi- 
mizing the  expected  distance  to  the  patrolling  force) .   In 
this  case  the  problem  takes  on  a  game  theoretic  aspect,  and 
results  from  game  theory  might  be  applicable. 

If,  through  careful  analysis  of  an  adversary's  intelli- 
gence gathering  capabilities,  force  levels  and  historic 
political  intentions,  one  is  able  to  determine  an  a  priori 
distribution  of  the  location  of  T,  then  it  is  possible  to 
use  the  probability  integral  transform  to  transform  the  a 
priori  distribution  into  a  Uniform  (0,1)  distribution.   This 
might  enable  one  to  utilize  the  models  developed  herein 
and  then  transform  the  results  back  to  the  original  a 
priori  space. 

C.   ECONOMIC  COST  OF  VARYING  THE  RADIUS  OF  INFLUENCE 

There  are  economic  costs  involved  in  varying  the  radius 
of  influence  of  a  carrier.   Through  the  use  of  the  radius 
of  influence  model  developed  herein  and  basic  cost  considera- 
tions, an  analysis  of  the  cost  tradeoffs  involved  in  achieving 
a  particular  improvement  in  P(D_<d)  could  be  developed.   This 
development  could  include  economic  tradeoffs  in  terms  of 
various  combinations  of: 

1.  Increasing  a  carrier's  radius  of  influence 

2.  Increasing  the  transit  speed  of  the  carrier 
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3.  Increasing  the  number  of  carriers 

4.  Increasing  -intelligence  concerning  the  location 
of  T. 
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